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introduction 
Problem  Statement 

TMs  fellowship  a«e.es~ 

:::::o«^Sof:ua,i«eape^so^^^^ 

physicists.  The  seventy  of  the  impact  of  m  ^  ju„ature  Many  approaches  to 
7uS  women  has  been  described  exsenswe^  m  ^  ^ 

the  improvement  in  the  overall  ,^3  Hem  is  re^ch  in  improved  breast 

being  pursued  Among  >Pf  ‘^'T^atter  sLses  the  need  for  the  general 

imaging  techniques  for  ^  ler  ®  techniques  to  monitor  the  population  at  nsk. 

availability  of  safe  and  low-cos  persistent  mortality  rate  despite 

This  is  particularly  urgent  in  view  of  explained  by  the  fact  that  minonty, 

improvements  in  diagnosis  and  treatment  is  m  part  t^plai  J  4 

older  and  low-income  women  have  tool  for  the  detection  of 

conventional  mammography  is  qui  e  ®  cancer  is  frequently  complicated  by  the 

microcalcifications,  early  detection  of  b  formation,  calcium 

occumence  of  benign  "  J„“tllW  hypetplasia  and  are  normal 

and  cholesterol  precipitation,  ‘7^”  ^f  these  benign  lesions  most  commonly 

variations  of  breast  tissue  fourth,  glandular  tissue  and  these  lesions  generally 

occurs  in  the  fa«y  evolution  of  the  breast 

^r^I^stp^rS  of  ."coirect  dia^osis  of  early  mammary  carcinoma. 


Background 

The  accreditation  of  medical  physicists  “ 

Mammography  by.rfe“'°“7Th^^”rLM  College  of  Radiology  (ACR)«,  the 
Physicists  m  Medicine  (A  )  ,  of  California  (Dept,  of  Health  Services) 

Federal  Drug  Administration  (  )  :,gj^ents.  The  former  were  addressed 

has  general  education  as  _  of  Radiological  Sciences  at  the  University  of 

through  graduate  studies  in  t  e  professional  training  courses. 

wrCwork  whhorcTrlipeLtion)  association  with  the  Veterans  Affairs  Medical  Center 

in  La  Jolla,  CA.  ♦  ;ms,crine  has  evolved  over  many  years  to  be  X-ray 

The  standard  for  diagnostic  breast  of  dusters  of 

“lSioi;i^rfondt:»rn^nesl’nbord.^ 


often  dffioult  to  distinguish  front  benign  fibrocystic  change  ^d  this  leads  to  a  high  number 
of  false  positive  diagnoses  and  thus  unnecessary  biopsies.  ’  There  has  been  extensive 
re^rch  te  provide  Idjunctive  or  alternative  imaging  systems  to  mammography 
digital  x-ray  imaging,  magnetic  resonance  imaging  (MSI)  and  ultrasound  sonography.  In  the 

eariy  1980’s,  dedicated  ultrasound  scanners  were  evaluated  for  whole  bre^  ii^ng. 

Today  the  role  of  ultrasound  in  breast  imaging  is  pnmanly  liimted  to  evaluating  the 
ICosWon  of  larger  lesions  detected  by  mammography  in  terms  of  cyst, c  (bemgn)  versus 

sohd  (CT)  is  a  powetfiil  and  routme  in  t-ra^  nuclear 

medtcine  and  magnetic  resonance  imaging.  In  the  1970’s  and  eariy  1980  s,  several 
“festigltors  sought  with  fimited  success  to  app^  these  image  reconstructton  metho  s  to 

pulse-transmission  ultrasound  of  breast  tissue.  The  goal  ^ 

improve  differentiation  of  breast  cancer  from  numerous  bemgn  renditions  of  the  br^st  by 
rruring  several  acoustic  properties  of  the  lesion,  A  particularly  ,n,er«,mg  approach  ,0 
ultrasound  CT  is  diffiaction  tomography  which  has  received  s^stantia  t  “■'n 
treatment  and  holds  the  potential  to  account  for  the  inherent  diffracuon  in  sou 
propagation.  Despite  considerable  theoretical  development  of  diffraction  tomograp  y, 

laboratory  implementation  has  been  limited  and  has  met  with  variable  . 

The  work  which  constitutes  the  dissertation  research  portion  of  tlus  feUowship  is 
based  ol^nl  approach  to  ultrasound  diffraction  tomography  which  is  largely  rile  resute 
of  a  wllaboration  between  researchers  at  UCSD  and  Thermotrex  Technologies  Co^ration 
(TTC).^‘'“  This  approach  is  described  briefly  in  a  subsequent  secuon  (Methods). 


Purpose 

Certification  of  medical  physicists  in  Diagnostic 
Mammography  permits  them  to  assure  and  verify  the  proper  and  safe 
diagnostic  imaging  equipment  and  thus  protects  the  safety 

m^tenance  of  imaging  performance  standards  results  in  improved  and  more  accurate 

diagnostics  and  thus  serves  the  general  patient  population.  This  “  “ 

mammography  where  the  imaging  procedures  have  to  be  optimized  to  permt  mites 
detectiof  of  breast  cancers.  This  has  recently  been  recognized  by  the  establi^ent  of  mo 
stringent  standards  both  for  the  certification  of  medical  physicists  md  the  prformance  of 
msimninaranhv  centers  Thorough  training  in  diagnostic  imaging  physics  and 

radical  phylcists  thus  proves  to  be  a  ben^t  to  the  local  commumty  in 
“e  fkiUs  and  training  acquired  will  serve  the  general  population  a,  nsk,  as  wefi  as  older, 
“o“^r  women  wSo  presently  receive  inadequate  diagnostic  screemng  for  breast 

The  principal  purpose  in  the  development  of  the  ultrasound  CT  scanner  by  the 
UCSD/TTC  team  was  that  the  differences  described  above  between  this  new  approach  and 
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previous  efforts  to  construct  a  practical  ultrasound  CT  (USCT)  system  might  prove  to  be 
advantageous^  The  principal  advantages  were;  1)  Low-ftequency  continuous-wave 
insonification;  2)  large  number  of  projections  around  the  whole  object,  3)  high  spee 
sampUng  at  optimal  spatial  frequency;  4)  large  field  of  view.  A  laboratory  prototype, 
described  in  the  Methods  section,  suggested  that  the  approach  may  have  utility  as  an 
adjunctive  tool  to  mammography.  The  human-sized  device  which  was  developed  produ^s 
unique  images  of  a  very  large  eight-inch  field  of  view  m  a  coronal  format.  In  cases  where 
x-ray  mammogram  reveals  a  suspicious  lesion,  monitoring  with  USCT  could  provide  a  non- 

invasive  alternative  to  biopsy.  ...  • 

The  focus  of  the  dissertation  research  is  the  quantitative  descnption  of  the  imagi  g 

characteristics  of  this  new  ultrasound  computed  tomography  procedure.  Diffraction 
tomography  potentially  affords  the  characterization  of  certain  breast  tissue  properties 
such  as  speed  of  sound,  density,  attenuation  and  compressibility.  Before  a  quantitative 
interpretation  of  such  parameters  can  be  realized,  however,  the  system  needs  to  be 
calibrated  with  phantoms  that  estabUsh  the  linearity,  resolution  and  dyn^c  range  of  the 
reconstructed  images.  This  is  particularly  essential  in  the  development  of  a  climcal  breast 
imaging  system  since  breast  tissue  is  complex,  the  breast  has  few  recognizable  anatomical 
landmarks  and  thus  imaging  artifacts,  noise  and  aberrations  must  be  carefolly 
characterized  and,  if  possible,  eliminated,  in  order  to  reduce  the  probabihty  of  mis¬ 
diagnosis  or  misinterpretation  of  the  images. 


Methods 

The  emphasis  of  this  pre-doctoral  fellowship  is  on  two  areas  that  are  relat^  in 
subject  matter  but  rather  distinct  in  terms  of  the  activities  and  approaches  involved^  The 
first  task  addresses  preparation  for  board  certification  in  Diagnostic  Radiological  P  ysics 
and  Mammography  and  thus  contains  both  formal  and  on-the-job  traimng  e  ements. 
Certification  in  diagnostic  radiological  physics  requires  both  formal  trainmg  leading  to  a 
graduate  degree  in  physics,  medical  physics  or  a  related  subject  area  as  weU  as 
documented  work-experience  in  the  field  under  the  supervision  of  certified  medical 
nhvsicists.  This  report  describes  the  training,  work  expenence  and  degree(s)  conterred 
Lring  the  first  year  of  this  fellowship.  The  second  task  involves  dissertation  research  on 
a  full-breast  imaging  technique  based  on  ultrasound  diffraction  tomography  pnnciples 
realized  in  a  recently  developed  scanning  system. 


1.  Certification/Diagnostic  Radiological  Physics  and  Mammography 

The  AAPM  Report  on  the  Clinical  Training  of  Radiological  Physicists  defines  the 
areas  of  activity  of  the  Diagnostic  Radiological  Physicist  as  follows;  1)  Calibration  of 

imaging  equipment;  2)  Calculation  and  measurement  of  exposure  and  dose;  3) 

Improving  and  maintaining  medical  image  quality;  4)  Training  of  physicists,  climcal 
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diagnostic  imaging  residents,  radiological  and  ultrasound  technologists,  and  other  allied 
heih  professionals  in  diagnostic  radiology;  5)  Education  of  health  professionals  and  the 
public  in  diagnostic  imaging  physics  and  radiation  effects;  6)  Additional  duties  such  as 

developmental  studies,  administrative  duties  and  research.  ^  .  ,  , 

Competence  in  the  above  radiological  and  clinical  areas  denves  from  fundamental 
understanding  of  and  facility  with  radiation  physics,  general  and  specific  imapng 
concepts,  radiation  biology  and  basic  clinical  concepts  including  anatomy,  physiology  and 

patient  American  College  of  Radiology  (ACR)  has 

established  qualification  requirements  for  accreditation  of  medical  physicists  in 
mammography  under  the  Mammography  Quality  Standards  Act  (MQSA).  These 

specifr  that  for  eligibility,  one  oftwo  sets  ofcriteria  must  be  met; 

1)  Medical  physicist  certified  in  Diagnostic  Radiological  Physics  or  Radiological 
Physics  by  any  certifying  body  recognized  by  the  American  Collep  of  Radiology. 

2)  Until  January  1,  1996,  alternative  criteria  will  be  recognized. 

a)  M.S.  or  M.A.  degree  in  physics,  applied  physics,  radiological  physics, 

biophysics,  health  physics  engineering,  public  health  when  B.S.  is  m  a  physical  science. 

b)  Training  in  biological  sdences. 

c)  Minimum  1  year  training  in  medical  physics  in  the  area  of  diagnostic 

radiological  physics.  l  ♦ 

d)  Minimum  2  years  experience  conducting  mammography  equipment 

performance  evaluations.  ^  .  u 

Additionally,  it  is  specified  that  a  candidate  satisfying  either  of  the  above 

requirements  have  had  at  least  fifteen  (15)  hours  of  documented  continuing  medical 
education  training  in  mammography  physics  within  the  last  three  years. 


2.  Dissertation  Research 

The  focus  of  the  doctoral  dissertation  is  the  quantitative  description  of  the 
imaging  characteristics  of  a  new  ultrasound  computed  tomography  procedure  that  is 
intended  to  serve  as  an  adjunct  to  conventional  x-ray  mammography.  This  new  approach 
is  based  on  two  human-size  systems  designed  specifically  for  breast  imagmg  that  use  low- 
frequency  (0.2- 1.4  MHz)  continuous-wave  sound  ,  cylindrical  arrays  of  1024  and  512 
transducers,  high  spatial  frequency  sampling  of  the  wavefront  and  a  unique  diffraction 
tomographic  reconstruction  method. 

This  method  of  ultrasound  CT  employs  a  fuU  circle  of  transducers  which 
surrounds  and  defines  the  field  of  view  (Figure  1).  The  object  is  coupled  to  the  transducer 
array  by  a  water  bath.  One  element  at  a  time  transmits  a  continuous  smgle-frequency 
sound  beam  with  a  cylindrical  wavefront.  The  system  reaches  steady  state  in  about  300 
psec,  at  which  point  the  remaining  transducers  record  amplitude  and  phase  (relative  to  the 
transmitter)  of  the  signal  scattered  the  object  to  be  imaged  .  Each  transducer  is  activated 
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as  a  transmitter  in  turn,  sequentiaUy  probing  along  a  complete  360 

The  laboratory  systems  were  constructed  to  test  feasibility  o  t  ®  ^  ° 
breast  imaging.  The  1  MHz  SCT  device  shown  in  Fi^re  2  consists  of  1024  lead 
zirconate  titanate  (PZT)  transducers,  arranged  in  an  eight-mch  cylinder^  The  breas 
imaging  sysiem.  shown  in  Figure  3,  consists  of  a  comrol  console  which  contains  all  the 
signal  generation  and  data  acquisition  and  processing  electromcs  and  computers  an 
couch  on  which  the  patient  reclines  with  the  breast  suspended  in  a  water  tank 
incorporating  the  transducer  array.  The  tank  can  be  raised  and  lowered  for  slice  selection 

and  moved  laterally  for  optimum  patient  positiomng.  ^ 

The  emphasis  of  this  year’s  work  was  on  evaluating  the  acoustiMl  charartenst  cs 

of  the  two  transducer  rings,  on  establishing  initial  system  petformant*  figures  “f 
on  assessing  the  effects  of  various  object  and  data  acquisition  variables  on  image  quality 

and  on  participating  in  a  clinical  trials  series.  These  tasks  followed  the  pta  in  the 

Statement  of  Lrk  of  the  original  proposal.  An  additional  task  was  mtroduc^.  a^lysts 
ofl  raw  sign^s  scattered  by  breast  tissue.  Tins  work  is  bas^  on  " 

analyses  reported  by  Steinberg  and  co-workers  as  well  as  other  investigators  who 
have  studied  the  amplitude  and  phase  distortions  of  the  wavefront  of  the  insomfyi,^ 
beam  when  passing  through  tissue  and  frequently  breast  tissue.  The  re  evance  o 
of  data  analysis  is  that  it  describes  and  summarizes  the  strength  of  the  t  e 

ultrasound  Lm  and  thus  was  felt  to  be  potentially  an  indicator  of  the  suitability  or  lack 
thereof  of  certain  reconstruction  algorithms  based  on  linearizing  assumptions  such  as  t 

^The  next  section  contains  a  summary  of  the  work  performed  during  the  year  and 
examples  of  some  of  the  results  obtained,  as  well  as  brief  descriptions  of  the  expenmental 
and  analysis  techniques  employed.  The  reporting  is  in  summary  form  since  most  results 
are  intermediate  and  incomplete.  The  tasks  addressed  are  summarized  below. 

.  Characterization  of  the  transmit  beam  patterns  and  slice  profiles  of  both 
the  1  MHz  and  0.5  MHz  SCT  rings. 


•  MTF  (Modulation  Transfer  Function)  measurements  vs.  frequency  and 
position  within  the  measurement  tank  for  the  1  MHz  SCT  ring. 

•  Speed  of  Sound  contrast  linearity  for  various  size  objects  in  a  non- 
aberrating  medium  (water)  for  both  1  MHz  and  0.5  MHz  SCT  rings  determined  wtth 
several  reconstruction  algorithms. 


Scattering  amplitude  distribution  studies  of  breast  tissue  from  four 
patients  scanned  during  the  clinical  trials  phase. 


Presentation  of  papers  at  the  AAPM  Annual  Meeting  in  Anaheim,  CA  in 
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July  1994,  at  the  American  Institute  of  Ultrasound  in  Medicine  Conference  m  San 
Francisco,’cA  in  March,  1995,  and  preparation  of  papers  to  be  presented  at  the  22nd 
International  Acoustical  Imaging  Conference,  Florence,  Italy,  Sept.  1995. 
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Figure  1.  Scanning  Geometry  of  Ultrasound  Diffraction  Tomography  Breast 
Imaging  System. 


•  1024  transducers,  20  cm  bore 

•  Operating  frequency, 700  kHz  -  1.2  MHz 

•  0.37  mm  image  resolution 

•  Data  acquisition  time,  0.5  sec 

•  image  reconstruction  time,  2  minutes 


Figure  2.  1024-element  1  MHz  SCT  Transducer  Ring 


CLINICAL  TEST  UNIT 


LO^ 


Figure  3.  TTC/LORAD  SCT  Clinical  Test  Unit 
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BODY 


This  section  describes,  largely  in  summary  form,  the  accomplishments  of  the  first  year  of 
the  pre-doctoral  feUowship  in  both  the  training  and  dissertation  research  areas.  M^y  results  of 
the  research  on  the  performance  evaluations  of  an  ultrasound  diffraction  tomography  breast 
imaging  system  are  preliminary  or  incomplete;  for  others,  the  requisite  data  analyses  have  not 
beerfi^ized.  The  thrust  of  this  report  then  is  not  to  present  final  results  and  conclusions  about 
the  technology  issues  raised  in  the  proposal  for  this  work,  but  to  indicate  progress  made  and  to 
point  out  areas  which  are  new  or  which  differ  significantly  from  the  statement  of  work. 

L  Certification  in  Diagnostic  Radiological  Physics,  Mammography 


Below  are  summarized  the  training  activities,  formal  education  elements  and  practical 
work  experience  that  all  satisfy  various  requirements  for  board  certification  in  Diagnostic 
Radiological  Physics  and  in  Mammography  by  agencies  such  as  the  American  Msociation  of 
Physicists  in  Medicine(AAPM),  the  American  CoUege  of  Radiology  (ACR),  the  Amencan 
College  of  Medical  Physics  (ACMP),  the  State  of  California  (Dept,  of  Health  Services)  and  the 

Federal  Drug  Administration  (FDA). 


1.1  Formal  Education 

Course  Work— Three  academic  quarters  of  Clinical  Rotation  in  several  areas  addressing 
various  imaging  modalities  were  completed  in  the  Radiological  Sciences  Department  of  the 
University  of  California,  Irvine,  CA  (UCI).  These  courses  were  conducted  mostly  at  the  UCI 
Medical  Center  in  Orange,  CA  but  also  included  visits  to  imaging  centers  such  as  the  UCI  Bram 
Imaging  Center.  These  rotations  emphasized  observation  of  and  participation  in  daily,  routine 
activities  in  a  hospital  setting  that  related  to  the  operation,  quality  control  and  interpretation  of 
various  imaging  systems.  The  rotations  were  as  follows: 

.  X-Ray  Radiography/Fluoroscopy-Participation  in  system  quality  control 

activities,  acquisition  of  images,  review  of  x-rays  with  radiologists. 

.  X-Ray  CT-Calibration  of  a  CT  Scanner  using  a  Calphan  phantom,  measurements 
of  system  resolution,  contrast  sensitivity,  slice  thickness,  linearity.  Observation  of  patient  scans 
and  review  of  tomographic  images  with  radiologists. 

•  MRI— Clinical  applications  presentation  by  radiologists,  observation  of  patient 
scans,  review  of  patient  images  with  radiologists,  daily  and  weekly  system  calibration  procedures 
using  several  phantoms. 

•  Mammography/Ultrasound-Participation  with  sonographers  on  patient  scans, 
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evaluation  of  mammograms  and  adjunctive  ultrasound  scans  of  a  range  of  breast  ^es  and 
pathologies  with  radiologists,  daily  ultrasound  and  mammography  systems  checks  an 

calibrations. 


.  Radiation  Oncology-Review  if  linac  calibration  procedures,  generation  of 

patient  depth-dose  profiles  from  bremsstrahlung  and  electron  beam  irradiation,  attendance  of 
weekly  tumor  review  boards,  review  of  treatment  planning  procedures  and  participation  in 

patient  treatments. 


.  Nuclear  Medicine-Presentation  of  and  participation  in  radiopharmaceuticals 

preparation  and  assaying,  observation  of  gamma  camera  patient  imaging,  daily  camera  calibration 
procedures,  review  with  radiologists  of  SPECT  images,  review  of  PET  procedures  mcluding 
cyclotron  operation  for  radiopharmaceutical  production,  observation  of  patient  treatment  and 
imaging,  review  of  images  with  radiologists. 


.  Digital  Angiography-Review  of  digital  fluoroscopy  procedures,  participation  m 

system  performance  calibrations,  implementation  of  various  digital  subtraction  techmques. 


Degrees  Conferred:  M.S.  degree  in  Radiological  Sciences,  March,  1995.  Umversity 
of  California,  Irvine,  CA. 


1.2  Training/Clinical  Experience 


.  Mammography  Practicum-Participated  in  a  1  week,  40  hour  practicum  for 

mammography  technologists  offered  by  the  UCSD  School  of  Medicine  at  the  UCSD  Medical 
Center  The  course  consisted  of  both  lectures  and  working  sessions  and  emphasized  breast 
pathology,  optimum  mammography  techmques  and  review  and  inteipretation  of  mammograms 

with  radiologists. 


•  Clinical  Trials— Participated  in  the  performance  of  clinical  trials  of  patients  at 
the  UCSD  Center  for  Women's  Health  designed  to  correlate  conventional  mammograms  and 
sonography  with  ultrasound  CT  images.  About  twenty  patients  with  a  wide  range  of  breast 
types  and  pathologies  were  studied. 


•  Diagnostic  Imaging/Mammography  Quality  Control-Participated  in  a  wde- 

ranging  number  of  diagnostic  imaging  equipment  quality  control  surveys  and  radiation  shielding 
planning  and  evaluation  procedures.  These  activities  are  summarized  in  a  table  in  the  Appendix^ 
Of  particular  importance  was  participation  in  the  acceptance  testing  of  a  GE  Seno^aphe 
Mammography  unit  at  the  Veterans  Affairs  Medical  Center,  La  Jolla,  CA.  in  Dec.  1993  and  the 
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subsequent  submission  of  calibration  and  performance  data  to  the  ACR  for  facility 
accreditation.  Accredidation  was  conferred  on  the  first  attempt  in  February,  1995. 


•  Mammography  Courses/Seminars— Attended  courses  dealing  with 
mammography  quality  control  practices  and  regulations  both  at  the  AAPM  Aimual  Meeting  in 
Anaheim,  CA  in  July,  1994  and  at  Southern  California  Chapter  AAPM  meeting  in  Carlsbad,  CA, 

May,  1995. 


•  AAPM  Summer  School-Received  29  CEU  units  in  Brachytherapy  Physics  by 
attending  the  1994  AAPM  Summer  School  at  UCSD,  San  Diego,  CA  in  June,  1994. 


H  DISSERTATION  RESEARCH 

Five  tasks  consistent  with  the  Statement  of  Work  for  Year  1  of  the  fellowship  proposal 
are  briefly  summarized.  Example  results  are  presented  and  some  preliminary  conclusions  are 
drawn  where  appropriate.  One  task  which  was  not  called  out  directly  in  the  proposal  is  the 
statistical  evaluation  of  the  raw  data  scattered  by  breast  tissue.  This  follows  recent  work 
reported  in  the  literature  where  such  analyses  are  used  to  correct  the  transducer  wavefi-ont 
distortions  in  conventional  B-scan  ultrasound  imaging.  The  application  here,  which  extends  the 
data  to  much  lower  fi'equencies  than  those  reported  before,  is  to  use  this  analysis  to  characterize 
breast  tissue  and  to  obtain  an  q  ptiori  indicator  of  the  degree  of  scattering  and  thus  the  validity  of 
particular  reconstruction  algorithms. 

During  the  reporting  period  a  second  transducer  ring  centered  at  0.5  MHz  was  fabricated 
and  integrated  with  the  system.  The  motivation  for  adding  such  a  low  frequency  ring  arose  out  of 
the  clinical  trials  which  indicated  substantial  attenuation  and  phase  shifts  for  either  large  or  very 
dense  breasts.  A  lower  fi-equency  ring,  though  resulting  in  loss  of  resolution,  was  expected  to 
have  better  signal  to  noise  characteristics  and  thus  provide  improved  images.  Many  of  the  tasks 
listed  in  the  Statement  of  Work  have  been  performed  for  both  transducer  rings,  thus  expanding 
the  level  of  effort. 

2.1  Transducer  Beam  Pattern  Characterization 

One  of  the  most  fundamental  system  characterizations  is  the  determination  of  the 
transducer  radiation  pattern  and  the  effective  slice  width.  The  radiation  pattern  must  be  specified 
accurately  for  solution  to  the  wave  equation.  For  the  sake  of  speed  of  calculation  and 
mathematical  simplicity,  it  is  presently  assumed  to  be  represented  by  the  first  derivative  of  the 
Hankel  function,  that  is,  it  is  assumed  to  be  dipole  radiation  with  a  cos0  distribution.  The  extent 
to  which  the  actual  radiation  (and  receive )  pattern  deviates  fi'om  such  a  distribution  is  of 
fundamental  interest. 

Measurements  were  carried  out  with  the  experimental  geometry  sketched  in  Figure  4.  A 
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“le  beam  profile  data  to  constant  radial  distance  from  the  transducer.  ^"51^ 
ZSucer  ring  the  needle  hydrophone  was  found  to  be  too  msensmve  to  record  the  SCT  (some 
XmS  rmography)  tran^ncer  signals  and  a  large,  0,24  in.  dtameter  transducer  was 

""’’’“’'Rgure  5  is  a  polar  plot  of  the  radiation  pattern  for  one  of  the  transducers  of  the  1  MHz 

SCT  ring  It  is  compared  with  both  a  cose  and  a  (cose  distribution  and  ts  seen  to  fit 

somewhere  inbetween.  The  situation  for  the  0.5  MHz  ring  is  d^ict^  m 

tha,  even  with  “tte  d"Son  Is^mpron  on  the  image 

to  width  (L/D)  ratio  (due  to  greater  detector  width,  1 .2  mm  vs.  .  mm 
in  a  non-cylindrical  antenna  pattern. 


NEEDLE  HYDROPHONE 


TRANSDUCER 

RING 


Figure  4.  Schematic  Diagram  of  SCT  Transducer  Beam  Pattern  and  Shce  Profile 

Determination  using  Needle  Hydrophone 
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Figure  5 .  Polar  Plot  of  Transducer  Radiation 
Pattern.  1  MHz  SCT  Ring. 


Figure  6.  Radiation  Pattern  for  0.5 
MHz  SCT  Ring  Showing  Side  Lobes. 
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2.2  Tomographic  Slice  Geometry 

The  measurement  of  effective  slice  geometry  was  performed  at  the  center  of  the  tank 
using  the  needle  hydrophone  for  the  1  MHz  ring  and  the  0.25  inch  dia.  transducer  for  the  lower 
frequency  ring.  One  transducer  was  excited  and  the  axial  pulse  intensity  distnbution  was 
measured  by  vertically 
translating  the  recording 
transducer.  The  results  for  the 
high  frequency  ring  are 
displayed  in  Figure  7  and 
suggest  that  the  slice  width  is 
essentially  that  of  the  transducer 
height  (L), 

though  the  profile  for  the  0.5 
MHz  ring  is  not  nearly  as  well 
defined  as  that  for  the  higher 
frequency  ring.  It  is  to  be  noted 
that  these  data  obtain  only  for 
one  transducer  at  one  location, 
the  center  of  the  tank. 


1.2 
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Figure  7.  Slice  Profile  for  1  MHz  SCT  Ring. 
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What  is  also  of  interest  is  the  effective  slice  width  or  sensitivity  profile,  measured  as  the 
axial  range  over  which  information  shows  up  in  the  reconstructed  image.  Such  a  measurement  was 
performed  with  two  orthogonal  nylon  filaments  inclined  at  approximately  45“  angles.  The 
filaments  were  120  pm  in  diameter  and  were  too  small  to  be  imaged  with  the  0.5  MHz  SCT  ring. 
With  the  1  MHz  ring,  however,  the  resultant  images  were  1.1  and  1.2  cm,  respectively,  very 
much  in  agreement  Avith  the  width  of  the  profile  in  Figure  7. 

2.3  System  MTF  Measurements 


A  phantom  incorporating  120  pm  diameter  nylon  filaments  arranged  in  rows  with  spacmg 
between  filaments  ranging  firom  0.5  mm  to  16  mm  was  constructed  and  used  to  measure  the 
point-spread-  function  (PSF)  at  various  locations  within  the  measurement  tank  as  well  as  at  a 
number  of  discrete  frequencies  within  the  envelope  of  frequencies  used  for  wide-band  data 
acquisitions.  These  measurements  were  carried  out  with  the  1  MHz  SCT  nng  m  water  and  thus 
represent  the  idealized  system  performance. 

In  general,  if  the  2-dimensional  system  point-spread  function  (the  image  of  a  point  object  or 
delta-function  input)  is  given  by  PSF(x,y),  then  an  object  with  true  spatial  characteristics  o(x,y) 
will  be  blurred  so  that  the  image  is  represented  by  the  convolution  integral: 


i{x,y)  =  JjJ dx'  dy'o{x'  ,y'  )PSFix  -x',y-y') 

which,  by  the  convolution  transform  theorem  is  equivalent  to  the  statement. 

l(,u,v)  =  CKu,v)*  H(u,v) 


(1) 


(2) 


where  I,  O  and  H  are  Fourier  transforms  of  i,  o  and  PSF  in  (1)  with  u,v  the  spatial  frequency 
variables  corresponding  to  x,y.  The  MTF  is  defined  by  the  ratio  of  the  blurred  to  the  true  object 

which  is  given  by: 

|6(M,t;)| 


MTF  =  - r 


=  H(u,v) 


(3) 


The  MTF  was  computed  from  measured  point-spread  functions  due  to  the  nylon  filaments  at 
several  locations  as  well  as  for  ten  discrete  frequencies  ranging  from  0.7  to  1.3  MHz.  A  typical 
result  of  such  measurements  is  shown  in  Figure  8a  which  suggests  that  the  cut-off  frequency  of 
the  MTF  is  somewhere  around  1  Ip/mm  which  is  roughly  on  the  order  of  the  theoretical  limit  for 
resolution  of  X/4.  We  did  not  find  much  variation  of  the  MTF  with  frequency,  however,  and  this 
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is  a  result  that  will  be  reviewed  further.  In  particular,  the  PSF  image  is  represented  by  veiy  few 
points  which  makes  the  MTF  calculation  (using  the  FFT)  less  accurate.  We  ^11  explore  the 
Lsibility  of  creating  a  speed  of  sound  phantom  with  a  sharply  defined  edge  for  whch  the^ 
spread-function  (LSF)  can  be  measured  at  some  angle  to  provide  more  data  points  for  the 

calculation,  a  method  commonly  employed  in  x-ray  CT. 

The  effect  on  the  system  resolution  of  the  presence  of  an  aberratmg  or  scattering  meium 
was  initially  explored  by  measuring  the  image  of  a  nylon  filament  passing  throu^  balloons  filled 
with  solutions  of  various  speeds  of  sound.  Also  shown  in  Figure  8b  is  the  MTF  for  such  a 
measurement  and  it  can  be  seen  that  the  cut-off  frequency  has  been  noticeably  reduced  from  that 
for  the  filament  in  water.  Thus,  it  is  expected  that  the  system  MTF  m  more  severely  scatteiing 
or  aberrating  media  such  as  agar  and  tissue-mimicking  phantoms  will  be  degraded  even  further. 
This  topic  will  be  revisited  in  experiments  planned  for  the  0.5  MHz  SCT  ring  as  well  in  the  near 

future. 


Figure  8  MTF  of  1  MHz  SCT  Ring.  First  Order  Bom  Magnitude  Reconstmction. 
Object  was  120  (tm  dia.  Nylon  Filament  (a)  in  Water,  (b)  inside  2.75  in.  dia.  Balloon  Filled  with 
Saline  Solution  with  Speed  of  Sound  4.7%  Higher  than  Water. 


2.4  Speed  of  Sound  Contrast  and  Linearity 

The  solutions  to  the  inverse  scattering  problem  based  on  the  Bom  approximation  result 

the  reconstmction  of  the  complex  propagation  constant. 
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(4) 


k\x,y)  =  ko(x,y{n\x,y)  + 1] 

where  ko  =  o)/Co;  n,  the  refractive  index,  is  given  by:  n(x,y)=c(x,y)/co  where  c(x,y)  is  the  speed 
of  sound  of  the  object  and  Co  that  of  the  medium.  Jhe  real  part  of  n  reconstructs  c(x,y),  the 
imaginary  part  is  an  attenuation  term  (i.e.  n(x,y)  -  a  +  ip  ). 

Measurements  of  the  linearity  of  the  reconstructed  speed  of  sound  over  a  range  "" 

-11%  to  +  9%  relative  to  that  of  water  were  earned  out  with  latex  balloon  phantoms  filled  wi 
various  ethanol  (lower  sound  speed)  and  saline  (higher  sound  speed)  solutions.  The  b^oons 
were  nominally  3/4  in.  in  diameter  thus  representmg  scatterers  of  dimensions  of  several 
wavelengths  at  both  nominal  SCT  ring  frequencies  of  0.5  and  1.0  MHz.  Simulation  calculations 
by  GP  Otto'  had  suggested  non-linear  behavior  of  the  reconstructed  sound  speed  and  attenuation 
les  for  ,he  variouf  re«,„s.ruCion  algorithms  currently  available.  The^  are  a 
algorithm  that  results  in  a  magnitude  image  (BFA)  and  two  approaches  to  comections  to  the  first- 
order  Bom  approximation  (SM  and  STT  algorithms)  the  discussion  of  which  is  mappropnate 

Measurements  were  conducted  with  both  SCT  rings  by  suspending  four 
and  obtaining  an  image.  For  each  image,  one  of  the  baUoons  was  filled  water  from  the  scanning 
tank  as  a  reference.  The  reconstructed  images  were  analyzed  by  recording  the  “"'ras*  (P‘“ 
values  relative  to  that  of  the  medium)  of  the  images  of  the  vanous  soluutms  for  both  the  real  md 
imaginaiy  images.  An  example  of  such  an  image  pair  is  shown  m  Figure  9  whch  depicts  tte  ted 
(9a)*lnd  imajiaiy  (9b)  parts  for  low-to-medium  comrast  solutions  and  qualitatively  satisfies  the 


Figure  9  Real  and  Itnaginaty  Components  of  Reconstmeted  Images  of  Medium  Contrast 
Balloons  using  Beam  Form  Algorithm  (TTC).  Clockwise  from  top  left:  Water,  1.73%  Saline, 

1.89%  SaUne,  -1.43%  Ethanol. 


'  GP  Otto,  Thermotrex  Technologies  Corporation,  San  Diego,  CA.,  Private  communication,  1995 
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satisfies  the  notion  that  higher  sound  speed  should  be  bright  and  that  the  attenuation  images 
should  all  appear  similar  since  no  significant  differences  in  the  attenuation  values  for  the  various 

solutions  were  expected. 


1  MHz  SCT  Ring-SOS  Calibration 
Reconstruction  Comparison-Real 


Figure  10.  Speed  of  Sound  and  Attenuation  Contrast  Sensitivity  for  1  MHz  SCT  Ring 
for  Three  Reconstruction  Algorithms.  Top:  Speed  of  Sound  (Real).  Bottom. 
Attenuation  (Imaginary). 
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Examples  of  speed  of  sound  and  attenuation  contrast  data  for  the  1  MHz  SCT  nng  are 
provided  in  Figure  10.  Similar  results  have  been  obtained  recently  with  the  0.5  MHz  nng  usmg 
the  same  sound  speed  solutions.  While  conclusions  are  premature  at  this  pomt  it  appears  that  as 
the  reconstruction  algorithms  become  more  sophisticated,  incorporatmg  time-of-flight  and 
multiple  scattering  corrections,  the  linearity  range  of  both  speed  of  sound  and  attenuation 
appears  to  increase  while  contrast  or  sensitivity  decreases.  It  should  be  noted  that  attenuation 
here  does  not  refer  to  actual  pressure  amplitude  or  intensity  reduction  but  simply  the  image  pixel 
values  for  the  imaginary  part  of  the  image  plotted  against  phantom  sound  spe^.  We  have 
already  initiated  experiments  to  measure  similar  data  when  the  speed  of  sound  objects  are 
embedded  in  strongly  aberrating  media  such  as  graphite-loaded  agar  and  foam  or  breast-tissue 
mimicking  materials  that  we  have  been  experimenting  with.  The  intent  here  of  course,  is  to  see 
how  the  speed  of  sound  contrast  and  linearity  are  affected  when  objects  of  dimensions  on  the 
order  of  several  wavelengths  are  embedded  in  materials  that  significantly  scatter  and  attenuate  the 

incident  wave. 


2.5  Scatter  Amplitude  Distribution  Analysis 


Wave  propagation  in  the  breast,  according  to  a  number  of  researchers,  is  ch^actenzed 

principally  by  scatter  and  refraction  and  these  processes  are  the  primary  causes  of  amplitude 
wavefront  distortion.  In  the  case  of  weak  scattering,  phase  distortion  can  be  adequately 
addressed  by  various  phase-aberration  correction  schemes.  Amplitude  distortion,  however,  is 
observed  when  scattering  becomes  strong  and  its  appearance  can  be  though  to  define  the  hmt  o 
applicability  of  weak-scattering  theories.  In  other  words,  scattenng  and  refraction  can  be  strong 
and  thus  one  would  not  expect  the  first  order  Bom  approximation  to  hold  smce  it  imphcitly 
assumes  that  scattering  is  weak.  The  evidence  for  this  transition  from  weak  to  strong  scattenng 
is  the  shape  of  the  observed  probability  density  function  (pdf)  of  the  complex  wavefront  the 
modulus  of  which  is  approximated  by  a  Rayleigh  distribution.  If  scattering  is  weak,  the  observed 

distribution  should  be  Rician  and  ultimately  Gaussian. 

Given  that  the  reference  waveform  in  the  homogeneous  medium  at  the  receiver  is. 

EJ^x)  =  A.{x)c\Tp{j(lf,{x)) 
and  that  with  the  object  in  the  path: 

E{x)  =  A{x)exp{j(f>(x)) 

the  wavefront  amplitude  distortion  term  is  defined  as: 
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(7) 


a{x)  =  A{x)l  \{x) 

Bv  formins  a  histogram  of  the  a(x)  values,  pdf  vs.  amplitude  plots  are  obtained  and 
comparS^to  Rayleigh  distLtions  of  identical  mean  values.  If  there  is  close  agreemen  between 
rZthatUd-the  standard  deviations  of  the  computed  Rayleigh  and  the  expenmentally 
»Ls  are  close,  then  the  scatter  field  U  considered  to  be  fully  developed. 

The  general  form  of  the  distribution  to  which  our  data  were  fitted  is  given  by  the 

expression:^^ 


A  \  A^  +  A^^J  ,^^o\ 
P(A)  =  ^exp,--^:;r-,4(  ^2  ) 


a 


r  2a^  J 


(8) 


where  P(A)  is  the  probability  density  for  amplitude  ^  ^  b^  intjM 

(8)  tends  toward  a  Rayleigh  distribution  whereas  for  large  values  of  A.,  P(A)  becomes  Oaussian. 

"  "™d“ri  of  an^ysis  to  lower  frequencies  by  studying  the  s^^^ 
amnUtude  distributions  from  four  patients  who  had  been  imaged  as  part  of  our  chmcd  mals 
series  at  the  UCSD  Center  for  Women’s  Health.  Two  patients  had  been  imaged  with  the  1  Mfc 
SCT  ring  and  two  with  the  new  low-frequency  ring.  One  sUce  from  each  patient  was  analyzed 
according  to  the  procedure  outlined  in  Appendix  2.  The  initial  analysis  used  only  those  reavers 
onhe  llsducer  ring  that  were  in  the  shadow  of  the  object.  This  is  essentiaUy  comparable  to  the 
techiique  used  by  Sttinberg  and  co-workers  who  used  a  single  transmitting  transducer  and 
linear  receiving  OTay  to  record  only  pulses  transmitted  through  the  object.  In  our  c^  we 
restricted  the  number  of  receiving  transducers  to  the  mimmum  value  detemuned 
examination  of  the  sinogram  of  the  data.  We  are  also  analyzmg  scatter  ampUtude  date  from  all 
receivers  surrounding  the  object  and  will  compare  that  with  the  sort  of  analysis  exemphfie  in 

ll  is  observed  that  the  index  is  lower  for  a  breast  that  had  had  a  lumpectomy  with 
attendanl  radiation  treatment  than  for  a  dei«  breast  (patients  5  and  8  -spectively)  and  to  the 
index  is  lower  for  a  fatty  breast  compared  to  a  normal  breast.  Since  these  data  are  from  only 
coronal  slice  it  is  not  clear  whether  this  ranking  will  hold  when  data  from  several  tomograp 
slices  of  each  breast  are  averaged.  While  the  behavior  with  frequency  for  the  low-frequency  rmg 
a^tLs  tot  cTl  to  cxp^tations,  the  classifying  index  based  on  the  mean  value  ^  standard 
deviation  of  the  recorded  amplitude  distribution  (i.e.  p  =  <x>/a  )  appears  to  differentiate 
between  breast  tissue  types. 
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BREAST  TISSUE  STATISTICS-SUMMARY 


FiKurell.  Scatter  Distribution  Index  u  =  <x>/a  as  a  Function  of  Frequency  for  Four 

Patients  with  Breast  Types:  Pat.  5;  Dense;  Pat.  8:  Fatty/Rad.Tx;  Pat.  14:NormaJ;  Pat.  15:Fatty 


2.6  Publications  and  Presentations 


Some  of  the  first  year's  work  of  this  fellowship  was  accepted  for  publication  or 
presentation  by  the  journals  and  conferences  listed  below.^^^^  This  effort  was  m  all  cases  the 
Llaboration  of  several  individuals  and  institutions  and  the  opportunity  to  participate  in  these 

efforts  is  gratefully  acknowledged. 


HS  Janee,  MP  Andre,  PJ  Martin,  GP  Otto,  BA  Spivey,  JP  Jones,  "Performance 
considerations  of  an  ultrasound  CT  system  for  breast  imaging".  Medical 
Physics  21(6);  1002,  1994. 

MP  Andre,  PJ  Martin,  HS  Janee,  GP  Otto,  BA  Spivey,  "Investigation  of  contrast 
sensitivity  and  phase  aberration  in  low-frequency  ultrasound  diffraction  tomograp  y  , 
Radiology  193(P);308,  1994. 
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JP  Jones,  HS  Jan^  “Quantitative  Comparison  of  the  Born  and  Rytov 
Approximations  in  the  Measurement  of  Scattered  Fields”,  to  be  presented  at  22nd  Int. 
Symp.  on  Acoust.  Imaging,  Florence,  Italy,  Sept.  1995. 

MP  Andre,  HS  Janee,  JP  Jones,  GP  Otto,  PJ  Martin,  “Reduction  of  Phase  Aberration  in 
a  Diffraction  Tomography  System  for  Breast  Imaging”,  to  be  presented  at  22nd  Int. 
Symp.  On  Acoust.  Imaging,  Florence,  Italy,  Sept.  1995. 

GP  Otto,  MP  Andre,  PJ  Martin,  HS  Janee,  BA  Spivey,  “Frequency-dependent  effects  on 
image  quaUty  in  a  diffraction  tomography  system  for  breast  imaging”,  Amer.  Inst. 
Ultrasound  in  Medicine,  San  Francisco,  March  1995. 


CONCLUSIONS 


Since  this  is  an  interim  report  that  covers  the  first  half  of  the  two-year  feUowship.  many 
of  the  conclusions  are  preliminaty  and  may  change  upon  further  mvesugattons  or  eMtnmattom  o 
present  data  and  resuhs.  Some  aspects  of  the  fellowship  grant,  m  parttcuta  those  deahng  wt 
training  and  certification,  require  no  conclusions  other  than  the  statement  that  certam  objecttves 
were  realized  and  others  are  left  to  be  accomplished. 


L  Certification  in  Diagnostic  Radiological  Physics/Mammography 


A  considerable  body  of  experience  in  diagnostic  imaging  equipment  quality  assurmce  and 
control  has  been  obtained  through  course  work  at  the  University  of  Califorma  Irvine  and  through 
surveys  of  many  systems  at  a  number  of  locations  under  the  supervision  of  a  board-certtfied 
medical  physicist.  Special  efforts  were  made  to  prepare  for  certification  in  mammography  with 
the  State  of  California  and  an  application  for  certification  wiU  be  submitted  in  the  near  foture. 

With  the  granting  oftheM.S.  degree  in  Radiological  Sciences,  expenence  gamed  under  the 

supervision  of  a  certified  physicist  can  be  counted  toward  the  work-expenence  requirements  for 
eUgibiUty  to  take  the  qualifying  board  examinations  offered  by  the  AAPM  or  ACR. 


BL  Dissertation  Research 


With  the  addition  of  the  low-frequency  transducer  ring,  many  of  the  performance 
ctaracterization  tasks  had  to  be  carried  out  for  both  units.  Thus,  the  Statement  of  Work  has  been 
expanded  in  scope  if  not  in  the  nature  of  the  tasks  proposed.  Ultimately  the  decision  may  be 
made  to  rely  on  only  of  the  transducer  rings  for  breast  imaging.  At  this  point,  not  enou^ 
information  is  available  about  the  relative  advantages  and  drawbacks  of  the  two  rings  to  formulate 

The  measurements  of  transducer  radiation  patterns  and  tomographic  slice  profiles  earned 
out  so  far  suggest  that  the  transducers  of  the  1  MHz  SCT  ring  produce  a  radiation  pattern  that  is 
a  fair  approximation  to  a  dipole  pattern.  This  cannot  be  said  for  the  lower  frequency  nng  the 
transducers  of  which  have  a  lower  L/D  ratio  and  thus  are  characterized  by  a  more  complex 
radiation  pattern.  It  will  be  one  of  the  tasks  for  next  year  to  investigate  the  effect  of  the  deviation 
from  a  dipole  pattern  on  the  reconstructed  images.  It  is  to  be  noted  here  that  the  measurements 
so  far  have  looked  only  at  a  few  individual  transducer  elements  and  that  it  is  assumed  that  the 
receive  pattern  of  the  transducers  is  similar  to  their  transmit  pattern.  There  are  indications  from 
the  examination  of  the  "noise  image",  that  is  the  image  of  the  water-filled  tank,  that  these 
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“'“"1re"Sckness  has  been  found  .o  be  close  to  the  long  dimension  of  the 

"^'em  MTF  measurements  to  date  have  been  only  earned  out  for  the  Wgher  frequency 

.ransducer  ring  and  suggestthat^^^ 

theoretically  expected  ”  feference’medium(water).  The  HTTP  becomes  degraded  when 

^obiS**sS^rsTStoder  filled  with  a  somewhat  different  sound-speed 

fhat  are  lot  yet  understood  but  may  reflect  a  basic  system  operation/image  reconstruction 
Speed  of  sound  contrast  and  linearity  measurements  have  been  ^ 

ijs,- . .. « j 

“EHSSESFzS:= 

Work  is  a  retsonable  outline  of  tasks  that  will  result  in  a  deeper  understanding  of  tte  quantitati 
performance  characteristics  of  a  new  breast  imaging  system  based 


scope  of  the  program  since  it  holds  the  promise  of  peimitung  the  optimiMhon  breast 
imln,  system  by  providing  frequency-dependent  tradeoffs  among  resolution,  contrast 
sensitivity  and  linearity  and  reduction  in  aberrations  and  image  artifacts. 
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appendix  1 


DIAGNOSTIC  IMAGING  QUALITY  CONTROL  EXPERIENCE  SUMMARY 


DATE 

location 

QCACnVITV 

8-7-92 

VAMC* 

GE  X-ray  CT  Unit  Annual  Survey 

9-8-93 

ff 

GE  Senographe  DMR  Mammography  Unit 
Acceptance  Testing 

6-23-94 

Children’s  Hospital** 

Mobile  X-Ray  Unit  Annual  Survey 

7-23-94 

Bi-Lateral  Cardiology  Unit  (Toshiba) 

Acceptance  Testing 

10-8-94 

ft 

Mobile  X-Ray  Units(5)  Annual  Survey 

12-4-94 

VAMC 

Accreditation  Testing,  GE  DMR 

Mammography  Unit 

12-15-94 

ft 

Cardiac  Cath  Lab  Radiation  Shielding 
Evaluation 

12-16-94  1 

Children’s  Hospital 

Radiography/Fluorscopy  Annual  Surveys, 

Rooms  2,3  and  5 

12-22-94 

ft 

X-Ray  CT  Scanner,  Room  4,5  Fluorscopy 

Units  Annual  Survey 

12-27-94 

VAMC 

GE  DMR  Mammography  unit  QC 

1-2-95 

Children’s  Hospital 

OR  C-Arm  Unit,  2  Mobile  C-Arm  units. 

Mobile  Dental  Unit  Annual  QC 

3-10-95 

VAMC  Outpatient 
Center*** 

Radiography/Fluoroscopy  Suite  Shielding 
Evaluation/Survey 

5-10-95 

VAMC 

Room  3  Fluoroscopy  Unit 

Annual  QC 

5-17-95 

ft 

Rpom  1  Radiography  Unit 

1  Annual  QC 
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*  Veterans  Affairs  Medical  Center,  La  Jolla,  CA. 


**  Sharpe/Children's  Hospital,  San  Diego,  CA 

***  Veterans  Affairs  Medical  Center  Outpatient  Clinic,  San  Diego,  CA 
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APPENDIX  2 


SCATTER  AMPLITUDE  DISTRIBUTION  ANALYSIS-BREAST  TISSUE 


The  histogram  data  of  the  ratio  of  scattered  pressure  amplitude  to  that  obsewed  with  the 
reference  medium  (water)  are  analyzed  to  see  whether  they  are  more  nearly  fit  with  a  Raylei^  or 
Rician  distribution  The  Rayleigh  distribution  in  scatter  amplitude  results  when  scattering  m 
reftaction  are  strong,  suggesting  that  the  object  is  composed  of  a  large  number  of  independent, 

random  scatterers  which  produce  unifoim  phase  distributions  and  noimal  amplitude 

distributions.  The  modulus  of  such  a  distribution  is  Rayleighian.  As  scattenng  and  refraction 
become  weaker,  the  distribution  tends  toward  a  Rician  distribution  and  ultimately  toward  a 

?rteuTindex  was  defined  as  :  p  =  <x>  /  o  where  <x>  is  the  mean  of  the  distribution  and 
o  its  standard  deviation  This  ratio  is  relatively  low  for  a  Rayleigh  distribution  and  increases  as 

the  distribution  tends  toward  Rician.  ....  ,  a 

An  example  of  a  simple  analysis  using  the  Mathsoft  Mathcad  apphcation  is  The 

example  is  for  a  tomographic  sUce  of  the  breast  of  Patient  5  which  was  classifi^  by  the 
mammography  radiologist  as  dense.  It  is  seen  that  the  scatter  amphtude  distnbution  best 
approximated  by  a  Rician  distribution  with  k  1 . 
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STATISTICAL  EVALUATION  OF  SCATTER  DATA  FROM 
BREAST  TISSUE 

PATIENT=5/SLICE  5  -  1.0  MHZ  SCT  RING  H.  S.  Janee 

FILE:  PAT53d4 

a.  Rflyli^igh  nistHhiition  Fit  to  Data 

m  =  0  ..  63 

=  READ  ( triads )  Read  in  Data  from  Triads  Files 

m  :=  32  ..63  max(x)  =  7.749*10^ 


-  Normalized  Histogram 

max(x) 

Determination  of  Std.  Deviation 


x1  =9.301 
x2  =  115.605 

o  =  Jx2  -  x1  ^  0  =  5.394 


63 


63 


^  (m-31)x^ 


^  (m  -  31 )  Xf, 


x1  = 


m  =  32 


63 

m  =  32 


X2 


m  =  32 


63 

Z 

m  =  32 


TISSUE  IMDEX: 


x1 

a 


1.724 


.  (m  -  31)- 


1 

■  ( m  -  31 )  f 

exp 

~2 

^ - —  1 

J  . 

3.254 
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